Robustness of Charge-Qubit Cluster States to Double Quantum Point Contact 
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We theoretically investigate the robustness of cluster states in charge qubit system based on 
quantum dot (QD) and double quantum point contact (DQPC). Trap state is modeled by an is- 
land structure in DQPC and represents a dynamical fluctuation. We found that the dynamical 
fluctuations affect the cluster states more than static fluctuation caused by QD size fluctuation. 
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One-way quantum computing [jj, y| is an important 
approach for quantum computation based on a series 
of one-qubit measurements starting from a cluster state 
of a qubit array. Cluster states are highly-entangled 
states involving all qubits and are typically generated 
from an Ising-like Hamiltonian, starting from an initial 
product state |*o) = |*(i = 0)) = n^i|+)j, where 
|±)i = (|0)j±|l)j)/\/2. Here, |0)i and |1), are the two 
states of the i-th qubit in an 7V-qubit system. In ref. Q , 
we showed that cluster states in charge qubits |4l-[l7| can 
be created by applying a single gate bias pulse, right after 
preparing the initial product state (one-step generation 
method), and are more robust against nonuniformities 
among qubits than decoherence-free (DF) states [1^ [l^] 
under a noise environment generated by a quantum point 
contact (QPC) detector, which is a sensitive detector of 
electric charge distribution [20, |21|. However, trap sites 
are often unavoidable in solid-state qubits owing to their 
small fabrication size and several experiments show that 
trap states significantly affect electric transport proper- 
ties of nanoscale devices |22h25| . Here, we model the trap 
site as the island (discrete energy state) between dou- 
ble QPCs (DQPC) and investigate robustness of cluster 
states in charge qubits measured by the DQPC detec- 
tor (Fig. 1). The charge qubits are based on quantum 
dots (QD), in which the position of the excess charge in 
a qubit affects the QPC current electrically, resulting in 
detection of charged state. Owing to this additional ex- 
ternal degrees of freedom, this setup is considered to be a 
harsher and more realistic condition for qubits than that 
of ref .[3]. Here, we calculate a time-dependent fidelity of 
four charge qubits with DQPC by solving density matrix 
(DM) equations. 



This article is organized as follows. In §2, formulation 
of our model is presented. Section 3 is devoted to the 
numerical calculations for the cluster state and DF state. 
The conclusion is given in §4 . The detailed coefficients 
used in the equations in the text are shown in Appendix. 




FIG. 1: Four qubits that use double QD charge states are 
capacitively coupled to a double QPC detector with an island. 
One excess charge is injected into each qubit. Single energy 
level is assumed in each QD and the island. This flgure shows 
1 1001} state. 



II. DENSITY MATRIX EQUATIONS 



Here we show the formulation of the density matrix 
equations for Fig. [T] The Hamiltonian for the com- 
bined qubits and the QPC for Fig. [T] is written as H = 
Hqh + Hqpc + Hint- ^qb IS & Hamiltoniau for an array 
of charge qubits with nearest-neighbor interactions, ex- 
pressed by 0, 0, i| 



-ffqb = ^(^l 0-ix + El (Ji 



7 , '^ij ^Jz "'jz ' 
i<j 



(1) 



where aix and aiz are Pauli matrices for the z-th qubit ex- 
pressed by a„ = |0),(1| + |1),(0| and a,, = |0),(0|-|1),(1|. 
|0) and |1) refer to the two single-qubit states in which 
the excess charge is localized in the upper and lower dot, 
respectively, fti is the inter-QD tunnel coupling between 
two coupled QDs. e^ is the charging energy and cor- 
responds to the energy difference between |0)i and |l)i 
controlled by gate voltage. The coupling constants Jy 
are derived from the capacitance couplings [26| . 

The Hamiltonian for the two serially coupled QPCs 



Hqpc is described by 



^qp<= =" X! X! [^^° ^L s C,„ s + Vi^ s (c|^ Js+dlc^^s) 






^ Eddld, + Ud\d^dld^ . 



(2) 



Here, q^s (ci^js) is the annihilation operator of an elec- 
tron in the z^th (iflth) level [iL{iR) — 1,2,3,...] of the 
left (right) electrode, ds is the electron annihilation oper- 
ator of the island between the QPCs, Ei^s {Eirs) is the 
energy level of electrons in the left (right) electrode, and 
Ed is that of the island. Here, we assume only one en- 
ergy level on the island between the two QPCs, with spin 
degeneracy. Vi^s iYiRs) is the tunneling strength of elec- 
trons between the left (right) electrode state ils {irs) 
and the island state. U is the on-site Coulomb energy of 
double occupancy in the island. 

TJint is the capacitive interaction between the qubits 
and the QPC, which induces dephasing between different 
eigenstates of aiz [la, [l^- Most importantly, it takes 
into account the fact that localized charge near the QPC 
increases the energy of the system electrostatically, thus 
affecting the tunnel coupling between the left and right 
electrodes: 
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where SVi^is (a = L,R) is an effective change of the 
tunneling strength between the electrodes and QPC is- 
land. Hereafter, we neglect the spin dependence of Vi^ 
and SVi^^i. We assume that the tunneling strength of 
electrons weakly depends on the energy Vi^i = Via{Ei^) 
and electrodes are degenerate up to the Fermi surface Ha ■ 
Then qubit states influence the QPC tunneling rate Pl 
and Tr by Tl^=T^^+T2^ and rj^^=r3"^-fr^^ through 

rf) = 2Trp^i^i^)\vt\t^»)\^ and PJ^^' = 27rp„(^„-H 
U)\V>^ {p,a+U)\'^, depending on the qubit state aiz = ±1 

i^L (M") = Viail^a) ± SVia{fJ.a) and pa{f^a) is the den- 
sity of states of the electrodes {a = L,R) and P^ is a 
tunneling rate when the island lies in the "c" state) . The 
values of P[ ^s are determined by the geometrical struc- 
ture of the system and depend on the distance between 
qubits and the QPC as P,^=^' = P,o± AP, = Po(l±C) with 
the measurement strength ( = AP/Pq |27| . 

The DM equations of four qubits and the DQPC de- 
tector at zero temperature of Fig. [1] are derived similarly 



to ref. [23, UM by 
dp^z. 



dt 

N 
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where z\ , zi indicate qubit states such as 
0000, 0001,..., 1111 and, p^^,,, p\z, and p\^z. 
are density matrix elements when no elec- 
tron ("a"), one electron ("b") and two elec- 
trons ("c") exist in the QPC island, respectively. 

•A)ooo = T^i e* + Jn + •^23 7 Jmm — X]j ^j ~ £4 + J\2 — J23: 
•■■, Jim = -J2i^i + di2 + J23- 9]{zi) is introduced 
for the sake of notational convenience and determined 
by the relative positions between qubit states (see 
Appendix and ref. [23|). In the present case, there are 
768 equations to be solved. 

Time-dependent fidelity F{t) = Tr[/5(0)/5(t)] is calcu- 
lated by tracing over the elements of the reduced DM 
obtained from eq. ^. 

We use these DM equations to describe four qubits, 
and compare fidelity of a cluster state 

I*)CS - ^(| + )l|0)2|+)3|0)4 + |+)l|0)2|->3|l)4 

+ h)l|l)2|-)3|0)4+h)l|l)2|+)3|l)4) (5) 

with that of a DP state 

l*>DF = ^(IllOO) - llOOl) - 10110) + 10011)). (6) 



III. CLUSTER STATES IN CHARGE QUBITS 

In one-way quantum computing, quantum computa- 
tion is carried out by measuring a series of qubits in clus- 
ter states. Cluster states are directly generated if we can 
prepare an Ising-like Hamiltonian Hcs = (ff/4)X]i<7(l ~ 
fiz)(l — o-jz), starting from the initial state |*o)- A uni- 
tary evolution Ucsit) = exp(— iiiJcs) (we use ?i = 1) at 
gt = (2n/ -|- 1)tt with an integer n/ transforms |^o) into 
cluster states. In the case of charge qubit in eq. ([1]), the 
aix term needs to be switched off during the creation 
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FIG. 2: Time- dependent fidelity of four-qubit states for a 
cluster state, a DF state, and a product state |1010), in case 
{!}. Fo = J, F; = Fi. The "weak" indicates a weak mea- 
surement case oi ( = 0.2 and the "strong" indicates a strong 
measurement case of ( — 0.6. 



of the cluster state and then switched on when mea- 
surements are carried out. To realize this controllabil- 
ity, additional gates are required. However, these addi- 
tional gates make the qubit system complicated, ending 
by producing decoherence and crosstalk between qubits 
themselves and between qubits and the environment. In 
addition, for some qubit systems, once the qubit array 
is made, fli and Jij are fixed, and only e^ is controllable 
via the gate voltage bias (we call these "simple-design 
qubits"). In ref. Q, we proposed how to effectively gen- 
erate cluster states in such simple-design qubits ("one- 
step generation method"). In this method, cluster states 
are obtained by applying a gate bias voltage Ci for the 
i-th qubit, expressed as 



Ef - nl 



(7) 



at time Us = vr(8n,7 + 1)/(4J). Here, Ei = (£» + 
ri^/ci) cosQfi with tanai — ^i/ei <C 1, and we also need 
a relation J (Sue — ni)/{8nj + 1) — Ei with an arbitrary 
integer nj and the number of the nearest qubits n^. As 
long as these equations are held, we can choose any inte- 
gers nj and fii. 

In the present one-dimensional qubit array, we have 
El = E^ = {She - 1)J and E2 = E3 = {Sue - 2)J. 
Here, we calculate two examples: case {1} is 17 = 2J 
{ue = 1) and case {2} is Q — iJ {he = 1). In case {!}, 



we have e^'' = £4'^ ~ 6.7J, e^^ — t 
we have e^*' — ef « 5.7J, 62^ — £3 



5.7 J. In case {2} 
4.5J. 



Figures[3] (a) and[3](b) show the combined effects of the 
island structure and the nonuniformities among qubits 
when parameters 17^, e; and F^ fluctuate by 10%. The 
smaller sizes of QDs are preferable for qubits based on 
QDsjTI, |8L26|. If we assume that a diameter of each QD 
is 5 nm[23], 10% size fluctuation corresponds to ±5 A size 
deviation from 5 nm. In both the strong measurement 
case [(a) C = 0.6] and the weak measurement case [(b) 
C = 0.2], the difference between the cluster state and 
the DF state is small. This is in contrast to the case 
without local island structure in QPC discussed in ref. [3|, 
and shows that the island structure (trap site) imposes 
a larger decoherence environment and would be a major 
origin of degradation for the entangled states. 

In order to examine the effect of nonuniformity be- 
tween qubits in more detail, we calculate a fldelity at 
time t — SOFq as a function of nonuniformity of qubit 
parameter 7/ (Fig. [4]). We can see that the difference is 
slight for various types of nonuniformities when rj is suf- 
ficiently small. We can also see that the fidelity of the 
larger gate bias case {1} is degraded shghtly more than 
that of the smaller gate bias case {2}. Note that this 
difference is much smaller in cluster states than in DF 
states, which is a preferable result for one-way comput- 
ing in solid-state qubits. From Figs. |3] and SI we can 
confirm that the cluster state is robust to nonuniformity, 
whereas the cluster state is weak against local electronic 
state (trap state), which is similar to the DF state. 
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FIG. 3: Time-dependent fidelities F(t) of a cluster state and a 
DF state when there are nonuniformities in qubit parameter, 
fi = 4J and ue = 1 (case {1}). Fo = J, F- = Fi, and C = 0.6. 
(a) Strong measurement case, (b) Weak measurement case. 
Nonuniformity in the qubit parameters is introduced as fii = 
2J(1— r^i), ei= e+riiJ, and F^ =(1— r;i)F' ', with i indicating 
the i-th qubit. Here rji = for all qubits besides r;4 = 0.1 (i), 
V2 = Vs = 0.1 (ii), and rj^ = 0.1 (iii). The fidelities of |*I')cs 
for (i) and (ii) mostly overlap. 



IV. NUMERICAL CALCULATIONS 



Figure [2] shows a time-dependent fidelity when the 
measurement strength is changed. It can be seen that 
when the measurement strength is sufficiently large (C — 
0.6), the fidelity of cluster states greatly degrades. This 
result shows that the cluster state is sensitive to the ex- 
istence of trap sites. 



V. CONCLUSIONS 

We investigated the effect of local electronic fluctua- 
tions (trap state) in addition to nonuniformities, as a 
model of a decoherence mechanism in cluster states. We 
found that the island (trap site) affects the fldelity of 
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FIG. 4: Fidelities of four-qubit cluster state and DF state at 
i = 50r^^ as a function of nonuniformity rj. Parameters are 
the same as those in Fig. [S] 



giizi)ii = 1, ..,4, zi = AA, ..., DD) are given by 



cluster states significantly. Because the local electronic 
state provides dynamical fluctuation whereas the nonuni- 
formities provide static fluctuation, we can say that the 
cluster states are sensitive to dynamical fluctuations. Ex- 
perimentally, it is difficult to reduce the number of trap 
sites. Thus, more detailed understanding of effects of 
trap states will be required in the future. 
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Appendix A: Coefficients in DM Equations 

Here, we show all coefficients used in eq. ^. For leg- 
ibility, we denote |00), |0f), |10), and |11) as \A), \B), 
\C), \D), respectively. Four-qubit states are written by 
\AA), \AB), ..,\DD). Then, J,^{zi = AA,...,DD) are 
expressed as 

JaA ~ £1+ £2+ £3+ £4+ J12-I- J23+ ^34, 

Jab ~ £1+ £2+ £3— £4+ J12+ ^23— J^34, 

J AC ~ £1+ £2— £3+ £4+ ^12+ ^23— ^34, 

J AD ~ £1+ £2— £3— £4+ J12— ^23+ ^34, 

JbA ~ £1— £2+ £3+ £4— J12— ^23+ ^34, 

JbB ~ £1— £2+ £3— £4— J12— ^23— ^34, 

JbC ~ £1— £2— £3+ £4— ^12+ ^23— ^34, 

JbD ~ £1— £2— £3— £4— J12+ ^23+ ^34, 

JCA ~ — £1+ £2+ £3+ £4— J12+ J23+ ^34, 

JCB ~ — £1+ £2+ £3— £4— J12+ J23— ^34, 

JCC ~ — £1+ £2— £3+ £4— J12— J23— ^34, 

JCD ~ — £1+ £2— £3— £4— J12— J23+ ^34, 

JdA ~ — £1— £2+ £3+ £4+ J12— J23+ ^34, 

Jdb ~ — £1— £2+ £3— £4+ J12— J23— J34, 

JdC ~ — £1— £2— £3+ £4+ J12+ J23— ^34, 

JdD = -£l- £2- £3— £4+ J12+ J23+ J34. (Al) 
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